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We present new measurements of time-dependent CP asymmetries for B° — > D(*)+I)(*± decays 
using (467±5) x 10 6 BB pairs collected with the BABAR detector located at the PEP-II B Factory at 
the Stanford Linear Accelerator Center. We determine the CP-odd fraction of the B° — ► D* + D*~ 
decays to be R ± = 0.158±0.028±0.006 and find CP asymmetry parameters S+ = -0.76±0.16±0.04 
and C+ = +0.00±0.12±0.02 for the CP-even component of this decay and S± = -1.80±0.70±0.16 
and C\ = +0.41 ±0.49 ±0.08 for the CP-odd component. We measure S = -0.63 ± 0.36 ± 0.05 and 
C = -0.07 ±0.23 ±0.03 for B° -> D+D~ , S = -0.62 ± 0.21 ± 0.03 and C = ±0.08 ± 0.17 ± 0.04 for 
B° -> D*+D~, and S = -0.73±0.23±0.05 and C = ±0.00±0.17±0.03 for B° -» D + D*~ . For the 
B° -> D*±pT decays, we also determine the CP-violating asymmetry A = ±0.008 ± 0.048 ± 0.013. 
In each case, the first uncertainty is statistical and the second is systematic. The measured values 
for the asymmetries are all consistent with the Standard Model. 
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I. INTRODUCTION 

In the Standard Model (SM), CP violation is described 
by the Cabibbo-Kobayashi-Maskawa (CKM) quark mix- 
ing matrix, V [1, 2]. In particular, an irreducible complex 
phase in the 3x3 mixing matrix is the source of all SM 
CP violation. Both the BABAR [3] and Belle [4] collab- 
orations have measured the CP parameter sin2/3, where 
(3 = a,Tg[-V cd V* b /V t dV* b }, in b -> (cc)s processes. 

The leading-order diagrams contributing to B° — > 
£)(*)+£)(*)- decays are shown in Fig. 1, where the color- 
favored tree-diagram of Fig 1(a) dominates. When ne- 
glecting the penguin (loop) amplitude in Fig. 1(b), the 
mixing -induced CP asymmetry of B° £)(*)+£)(*)-, de- 
noted S, is also determined by sin2/3 [5]. The effect of 
neglecting the penguin amplitude has been estimated in 
models based on factorization and heavy quark symme- 
try, and the corrections are expected to be a few percent 
[6, 7]. Large deviations of S in B° -> £)«+£)(*)- de- 
cays with respect to sin2/3 determined from b — > (cc)s 
transitions could indicate physics beyond the SM [8-10]. 

The CP asymmetries of B° — ► Z)W + decays have 
been studied by both the BABAR [11, 12] and Belle [13- 
15] collaborations. In the SM, the direct CP asymmetry 
C, defined in Sec. IV, for the B° -> £>(*)+£>(*)- decays 
is expected to be near zero. The Belle Collaboration has 
observed a 3.2 sigma deviation of C from zero in the 
B° — > D + D~ channel [15]. This has not been observed 
by BABAR nor has it been seen in other B° -> DW+DW- 
decay modes, which involve the same quark-level dia- 
grams. As was pointed out in [9], understanding any 
possible asymmetries in these decays is important to con- 
straining theoretical models. 




(a) Tree (b) Penguin 



FIG. 1: Leading-order Feynman graphs for the B { 
£>{*)+£>{*)- decays. 
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In this article, we update the previous measurements 
of CP asymmetry parameters in B° — > £>M+£)W~ de- 
cays [11, 12], including the CP-odd fraction for B° — > 
D* + D*~, using the final BABAR data sample. Charge 
conjugate decays are included implicitly in expressions 
throughout this article unless otherwise indicated. 



II. DETECTOR, DATA SAMPLE, AND 
RECONSTRUCTION 

A. The BABAR detector 

The data used in this analysis were collected with the 
BABAR detector [16] operating at the PEP-II B Fac- 
tory located at the Stanford Linear Accelerator Center 
(SLAC). The BABAR dataset comprises (467 ± 5) x 10 6 
BB pairs collected from 1999 to 2007 at the center-of- 
mass (CM) energy = 10.58 GeV, corresponding to 
the T(45) resonance. We use GEANT4-based [17] Monte 
Carlo (MC) simulation to study backgrounds and to val- 
idate the analysis procedures. 

The asymmetric energies of the PEP-II beams provide 
an ideal environment to study time-dependent CP phe- 
nomena in the B° — B° system by boosting the T(4S) 
in the laboratory frame, thus making possible precise 
determination of the decay vertices of the two B me- 
son daughters. BABAR employs a five-layer silicon vertex 
tracker (SVT) close to the interaction region to provide 
precise vertex measurements and to track low momentum 
charged particles. A drift chamber (DCH) provides excel- 
lent momentum measurement of charged particles. Parti- 
cle identification of kaons and pions is primarily derived 
from ionization losses in the SVT and DCH and from 
measurements of photons produced in the fused silica 
bars of a ring- imaging Cherenkov light detector (DIRC). 
A CsI(Tl) crystal-based electromagnetic calorimeter en- 
ables reconstruction of photons and identification of elec- 
trons. All of these systems operate within a 1.5 T su- 
perconducting solenoid, whose iron flux return is instru- 
mented to detect muons. 



B. Candidate reconstruction and selection 

The candidates used in this analysis are formed from 
oppositely charged mesons where we include the 

D*+ decay modes D*+ -> D°7r+ and D*+ -> D+ir° and 
D decay modes D° -> K~it + , D° -> R-tt+t: , D" -> 
K-TT+TT-TT+, D° -> K°tt+tt-, and D+ -> A'-7r+7r+. 
In the B° — > D* + D*~ mode, we reject B° candidates 
where both D* mesons decay to Dtt q because of its 
smaller branching fraction and larger backgrounds. Ref- 
erence [18] contains the details of the reconstruction pro- 
cedure, outlined here, used to select signal candidates. 
Charged kaon candidates must be identified as such using 
a likelihood technique based on the opening angle of the 
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Cherenkov light measured in the DIRC and the ioniza- 
tion energy loss measured in the SVT and DCH [16]. We 
reconstruct K® candidates from two oppositely charged 
tracks, geometrically constrained to a common vertex 
and with an invariant mass within 20 MeV of the nom- 
inal value [19]. We also require that the % 2 probability 
of the vertex fit of the be greater than 0.1%. We 
form 7T° candidates from a pair of photons detected in 
the calorimeter, each with energy greater than 40 MeV. 
The invariant mass of the two photons must be less than 
30 MeV/c 2 from the nominal 7r° mass, and their summed 
energy must be greater than 200 MeV. In addition, we 
apply a mass constraint to the tt° candidates. We re- 
quire the reconstructed D meson candidate mass to be 
within 20MeV/c 2 of the nominal value, except for the 
D a — > K~ir + ir° decays where we use a looser require- 
ment of 40 MeV/ c 2 . The daughters of each D candidate 
are fit to a common vertex with their combined mass con- 
strained to that of the D meson. We use D candidates 
combined with a pion track with momentum less than 
450 MeV/c in the CM frame to form D* + candidates. We 
fit the B° decay with a vertex constraint. 

Since the time of our previous publications [11, 12, 18], 
the BABAR reconstruction routines have been extensively 
revised, leading to significant improvements in localiz- 
ing and reconstructing tracks, particularly for low mo- 
mentum charged particles. These improvements have in- 
creased the reconstruction efficiency for final states with 
multiple slow particles, such as the B° — > D* + D*~ chan- 
nel which has a better than 20% improvement. As a 
result, the statistical sensitivity of the measurements in 
this paper has increased more than would be expected 
by just the increment in luminosity. 

To suppress e + e~ — > qq (q = u, d, s, and c) continuum 
background, we exploit the spherical shape of BB events 
by requiring the ratio of second to zeroth order Fox- 
Wolfram moments [20] to be less than 0.6. We select the 
B° candidates based on four variables: AE = E B — ^/~s/2, 
where E B is the energy of the B meson in the CM frame, 
the D candidate flight length significance, defined as the 
sum of the two D candidate flight lengths divided by the 
error on the sum, a Fisher discriminant [21], and a mass 
likelihood of the £>M mesons. The Fisher discriminant 
is a linear combination of 11 variables: the momentum 
flow in nine concentric cones around the thrust axis of the 
B° candidate, the angle between the thrust axis and the 
beam axis, and the angle between the line-of-flight of the 
B° candidate and the beam axis. The mass likelihood is 
formed from Gaussian functions, 



G(m D ;m DrDG ,a mD ) x G(m n ;m 



DpDG ' 



(Tm E ) 



/ core G(Ara D .+ ;Ara D .+ ,CAm 00 rJ 

PDG 

+ (1 - fcore)G(Am D *+; Am D .+ G ,(TAm te ii! 

x / CO roG(Ara£).-;Am D .- , <J Am core ) 

PDG 

+ (1 - /corc)G(Am £1 »- ; Am D - , 0-Amt«ii> 




(a) B° 



5.27 5.29 

m EI (GeV/c 2 ) 

D*+D*- 



(b) B° 



5.27 5.29 

m EI (GeV/c 2 ) 

D+D- 




5.23 5.25 

(c) B° -> D*+D 



FIG. 2: Projections of the tties fit results. The solid line 
represents the total fit PDF and the dashed line is the back- 
ground contribution. 



where the PDG subscript refers to the nominal value [22]. 
The reconstructed masses and uncertainties a m _ for the 
D mesons prior to the mass constraint are used in the 
likelihood. The D* portion of the likelihood is the sum 
of two Gaussian functions, a central core and a wider 
tail. The value of f COIC and the widths of the D* Gaus- 
sian functions are taken from detailed signal MC stud- 
ies, which show good agreement between data and MC 
samples. The selection criteria are optimized for each D 
decay channel to maximize the total signal significance 
S/ V S + B for each B° decay mode, where S and B are 
the signal and background yields, respectively. The opti- 
mized selections are specified in [18]. We keep candidates 
yfs/A-p*^ > 5.23 GeV/c 2 , where p* B is the 

On 



(1) 



With TOES 

momentum of the B candidate in the CM frame, 
average 1.1-1.8 candidates per event satisfy all of the se- 
lection criteria depending on the process. When more 
than one B° candidate meets the selection criteria, the 
one with the best C mass is kept. We find from MC that 
this procedure retains the correct candidate more than 
95% of the time. 

To determine the signal yields of the data sample, 
we use unbinned maximum likelihood (ML) fits to the 
toes distributions. The signal is described by a Gaus- 
sian function and the combinatorial background by a 
threshold function [23]. In detailed MC studies of the 
background, we find that there is a background contri- 
bution that exceeds the threshold function in the region 
toes > 5.27GeV/c 2 , where most of the signal events lie. 
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We describe this component with a Gaussian function 
having the same mean and width as the signal and re- 
fer to it as peaking background because if neglected, it 
would lead to an overestimate of the signal yields. In 
the B° — ► D* + D*~ channel, the peaking background 
arises primarily from misreconstructed B + — ► D* + D*° 
events where the slow 7r° from the D* a — > D°ir° de- 
cay is replaced by a n~ to form a D*~ candidate. For 
the other three processes, our studies of the composi- 
tion of the peaking background show it to be consis- 
tent with that of the combinatorial background in the 
region toes < 5.27GeV/c 2 . We treat the peaking back- 
ground component as an extension of the combinato- 
rial background. The peaking background yields rela- 
tive to the signal are fixed from MC to (1.6 ± 1.9)%, 
(7.1 ± 5.9)%, and (7.4 ± 2.9)% for the B° -> D*+D*-, 
B° -> D+D-, and B° -> D^D* modes, respectively, 
where the errors are due primarily to the size of the 
MC sample available for background studies. The sig- 
nal mean and background shape are free parameters in 
the fits. We fix the width of the signal Gaussian shape 
for B° -> D+D- and B° -> D^D* to 2.46 MeV/c 2 and 
2.55 MeV/c 2 , respectively, determined from MC, while 
the width of the B° — > D* + D*~ signal is allowed to float 
because of its much higher purity. The signal yields are 
934 ± 40 B Q -> D*+D*~ events, 152 ± 17 B° -> D+D- 
events, 365 ± 26 B° -> D*+D~ events, and 359 ± 26 
B° — > D + D*~ events, where the uncertainty is statis- 
tical only. The signal yields are consistent with previ- 
ously measured B° — > decay branching frac- 
tions from BABAR [18] and Belle [15, 24]. When com- 
pared with past BAB4R measurements [18, 25, 26], the 
low B° — ► D + D*~ yield in Ref. [11] is consistent with a 
statistical fluctuation. The fit projections for each mode 
onto Toes are shown in Fig. 2. 



III. TIME-INTEGRATED MEASUREMENT OF 
THE CP-ODD FRACTION 



The B° — > D* + D*~ process has two vector mesons in 
the final state and is an admixture of CP-even and CP- 
odd states depending on the orbital angular momentum 
of the decay products. We measure the CP-odd fraction 
R± using a time-integrated angular analysis [27] . We de- 
fine the three angles in the transversity basis as depicted 
in Fig. 3: the angle 9\ between the slow pion from the 
D*~ and the direction opposite to the D* + momentum in 
the D*~ rest frame; the polar angle 9 tr and the azimuthal 
angle tr of the slow pion from the D* + in the D* + rest 
frame where the z axis is normal to the D*~ decay plane 
and the x axis is opposite the D*~ momentum. Work- 
ing in the transversity basis, the time-dependent angular 




FIG. 3: Depiction of the B° ^_D* + D*~ decay in the 
transversity basis with the D*~ — > D°tv~ decay plane shown. 
The three transversity angles are defined in the text. 



distribution of the B° decay products is 
1 d 4 r 9 1 



T d cos 0i d cos 6> tr d^ tr di 16tt | A 1 2 

1 2 cos 2 9 1 sin 2 9 tI cos 2 (j) tI | A \ 2 

+ sin 2 6 1 sin 2 tr sin 2 4> tI \Au | 2 
+ sin 2 9 1 cos 2 6> tr |A ± | 2 
— sin 2 9i sin 20 tr sin0 tr ImfAjjAjJ 
1 



l^ll 



—= sin 29 1 sin 2 9 tr sin 2<p tt Re(A% A\\ ) 
v2 

—= sin 29i sin 20 tr cos (f>t T Ira(AQA±) 
v2 



(2) 



where Ak, with k = |,0, _L, represent time-dependent 
amplitudes given by 

_ V2A k (0) imt t/2To0 
Ak(t) - l + |A fe | 2e 

is — h tricp x k sin — ^— \ ■ (3) 

Here, t]q p is the CP eigenvalue, +1 for A^q, —1 for A±; 
Afe is the CP parameter defined in Sec. IV; Ato^ is the 
B° mixing frequency, (0.507± 0.005) ps _1 ; and t b o is the 
B° lifetime, (1.530 ± 0.009) ps [19]. Expressions similar 
to Eq. 2 hold for B° decays where each Ak is replaced by 
the appropriate A^ including A± — > — Aj_. Integrating 
Eq. 2 over t, </> tr , cos 9\ and averaging over B flavor while 
taking into account detector efficiency yields 



1 dr 



rdcos0 tr 32tt 

1 + a 



-^-(l-i? ± )sin 2 tr 



J (COS tl ) 



1 -a 



: J||(cOS0tr' 



-R_l cos 2 tr x 2x(cos0 tr ) : 



(4) 
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where we define 

R± = 

a = 



2 



l^jjl 2 

Kl 2 



\A° ± 



l^l 2 + KI 2 



and A® — Ak(0). The three efficiency moments 
/fc(cos^tr) are defined as 



4(cos( 



where go 



dcOS0id(£tr5fe(01 5 0tr)£(01,0tr, 0tr) , (5) 



4 cos 2 9i cos 2 



<Ptt, 9\\ 



2 sin 2 9i sin 2 tr , g± 



sin 0i , and e is the detector efficiency. The moments Ik 
are parameterized as second-order even polynomials in 
cos #tr whose parameters are determined from signal MC 
simulation and fixed in the fit. The three Ik functions de- 
viate only slightly from the same constant, making Eq. 4 
nearly insensitive to a, which we fix to zero in the fit. 

Because cos 9 tT is defined with respect to the slow pion 
from the D* + decay, the measurement resolution smears 
its distribution. We convolve the function from Eq. 4 
with a resolution function TZ(A9tr) which is modeled as 
the sum of three Gaussian functions. In addition, we 
include an uncorrelated Gaussian shape centered at 7r/2 
and normalized in < 9 tI < it to describe decays where 
the slow pion is poorly reconstructed leading to a loss of 
angular information. The uncorrelated term represents 
3% of the signal events where both slow pions are charged 
and around 16% in the modes where one of the slow pions 
is neutral. We determine the parameters of the resolu- 
tion model and of the uncorrelated term from signal MC 
simulation and fix them in the ML fit. Small differences 
observed in the angular distributions based on the charge 
of the slow pions lead us to divide the efficiency moment 
and resolution parameters into three categories, 7r°7r~, 
7r + 7T°, and 7r + 7r~. 

We determine R j_ in a simultaneous unbinncd ML fit to 
the mES and cos 9 tl - distributions for the three slow-pion 
modes. The iwes probability density function (PDF) was 
described in Sec. II B. The signal cos#tr distribution is 
given by Eq. 4 convolved with the resolution model. The 
background cos#tr distribution is modeled as a second- 
order even polynomial /b g (cos 9 tI ) = 1 + &2 cos 2 6*t r , where 
62, common to the three slow-pion modes, is allowed to 
float. The yield for each of the three slow-pion modes is 
determined by the fit. We validate the fitting procedure 
using high-statistics MC samples divided into data-sized 
subsets and find no significant bias. Fitting the data and 
including systematic uncertainties described below, we 
find 



R x = 0.158 ± 0.028(stat) ± 0.006(syst) . 



(6) 



Figure 4 shows the projection of the fit result. 

To evaluate the systematic uncertainty of R± , we vary 
the parameters used to model the efficiency moments 




FIG. 4: Projection of the fit result onto cos# tr for events with 
niES > 5.27GeV/c 2 . The solid line is the projected fit result. 
The dashed line is the background component. 



TABLE I: Summary of systematic uncertainties on the mea- 
surement of R±. 



Angular efficiency moments 
Angular measurement resolution 
a parameter uncertainty 
Peaking background 
cos# tr background shape 
Potential fit bias 



0.0024 
0.0036 
0.0026 
0.0014 
0.0002 
0.0017 



Total 



0.0055 



within the uncertainties of the MC simulation used to 
extract them. We do the same for the parameters used 
to model the experimental resolution. In both cases, 
we take into account correlations among the parameters 
when perturbing the values. We fix a to zero in the 
nominal fit, so we also set it to ±1 and assign the ef- 
fect on the fitted result as a systematic uncertainty. We 
change the mES and cos6* tr shapes of the peaking back- 
ground and assign the corresponding changes in R± as a 
systematic uncertainty. We allow the cos 9 tT background 
to have an additional fourth-order term to test our as- 
sumption of this background shape. This term is found 
to be consistent with zero, and we take the difference in 
R± with respect to the nominal second-order background 
description as the uncertainty with this model. We in- 
clude as a systematic uncertainty the statistical uncer- 
tainty associated with the MC validation. A summary 
of the systematic uncertainties is found in Table I. The 
total systematic uncertainty is the sum in quadrature of 
the individual contributions. 
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IV. TIME-DEPENDENT CP MEASUREMENT 

The decay rate /+(/_) of the neutral B meson to a 
common final state accompanied by a B° (B°) tag is 

f±(At) oc e HAt|/T *° {(1 T Aw) ± (1 - 2w)x 

[5 sin( Am d At) - C cos(Am d At)} } , (7) 

with CP asymmetry parameters S = 2Im(A)/(l + |A| 2 ), 
C =_(1 - |A| 2 )/(1 + |A| 2 ), and A = (q/p)(A/A), where 
A (A) is the decay amplitude for B° (B°) and q/p is 
the ratio of the flavor contributions to the mass eigen- 
states [28] . The parameter w is the average mistag proba- 
bility, and Aw is the difference between the mistag prob- 
abilities for B° and B° . Here, At = t ICCO — t tag is the 
proper time difference between the B reconstructed as 
B° -> DW+DW- (B rcc ) and the B used to tag the fla- 
vor (i? ta g). In the case of B° — > D* + D*~ , we obtain an 
expression similar to Eq. 7 from Eqs. 2 and 3, 

f±(At, cos(9 tr ) oc e - |At|/T B° {F(l t Aw) ± (1 - 2w)x 

[Gsin(Am d At) - H cos(Am d At)]} . (8) 

The F, G, and _ff coefficients [29] are 

F = (1 - R ± ) sin 2 6> tr + 2R ± cos 2 6> tr , 

G = (1 - R±)S + sin 2 6> tr - 2i? ± S'_ L cos 2 tr , 

= (1 - i?_L)G+ sin 2 <9 tr + 2i? ± G ± cos 2 6> tr . (9) 

The Afc parameters in Eq. 3 need not be the same be- 
cause of possible differences in the relative contribution 
of penguin and tree amplitudes, therefore the S and C 
parameters for each of the three (0, ||, _L) amplitudes can 
also differ. Note that the minus sign before S± in the 
expression for G absorbs rj^ p . We then define 

_ S ll \A^ + S \A^ _ q||Ajj| 2 + C |Ai]| 2 

|Ajj| 2 + |A[]| 2 ' G+ L4|j| 2 + |A[]| 2 ' 

(10) 

where A° k = A k (0) from Eq. 3. 

In the absence of penguin contributions, S D + D - = 
S+ = S ± = -sin2/3, and C D+D - = C+ = C ± = 
0. Because B° — > D* ± D T is not a CP eigenstate, 
the expressions for S and C are related, S d *± D t — 
—VI — Cd*± £it sin(2/3 c ff±(5), where 5 is the strong phase 
difference between B° -> D*+D~ and B° -» D+D*- de- 
cays [30]. Neglecting the penguin contributions, /3 e ff = /3, 
and C D »+ D - = -C D + D *-. 

The technique used to measure the time-dependent CP 
asymmetry is discussed in detail in Ref. [31]. We calcu- 
late At between the two B decays from the measured 
separation Az of their decay vertices along the z axis. 
The B rcc decay vertex is determined from the daughter 
tracks of the B° -> £>(*)+£)(*)- decay. The B tag decay 
vertex is determined in a fit of the charged tracks not be- 
longing to -B re c to a common vertex with a constraint on 
the beamspot location and the B rcc momentum. Events 



that do not satisfy |At| < 20 ps and cta* < 2.5 ps are 
considered untagged in the time-dependent fit. 

The flavor of the B tag meson is determined using a 
multivariate analysis of its decay products [31]. The 
tagging algorithm classifies the B flavor and assigns the 
candidate to one of six mutually exclusive tagging cat- 
egories based on the output. A seventh untagged cat- 
egory is for events where the flavor could not be de- 
termined. The performance of the tagging algorithm, 
its efficiency and mistag rates, is evaluated using the 
time-dependent evolution of a high-statistics data sam- 
ple of T(45) — ► -Btag-Bflav, where the -Bfl av meson de- 
cays to a flavor eigenstate D^~h + and h + may be a 
7T+, p + , or af . The tagging algorithm has an efficiency 
£tag = (74.4 ± 0.1)% and an effective tagging power 
Q = e t ag(l - w) 2 = (31.2 ± 0.3)%. The finite resolution 
of the B vertex reconstruction smears the distributions 
described in Eqs. 7 and 8. This measurement resolu- 
tion is modeled as the sum of three Gaussian functions 
described in Ref. [31], the parameters of which are also 
determined from the Bflav sample. 

We determine the CP asymmetry parameters in un- 
binned ML fits to the toes, At, and in the case of 
B° -» D*+D*-, cos6» tr distributions. The At signal 
distributions are given in Eqs. 7 and 8 convolved with 
the experimental resolution. The At background distri- 
bution has both zero and nonzero lifetime components 
which are convolved with the experimental resolution. 
The lifetime component is allowed to have effective CP 
parameters and lifetime, which are determined in the fits. 
The angular measurement resolution, determined for the 
CP-odd fraction measurement, is convolved with the sig- 
nal angular distribution. The efficiency moments are not 
modeled but rather absorbed into an effective R±, which 
is determined in the fit. This procedure simplifies the 
cos#tr distribution and does not introduce a bias. The 
peaking background for the B° — ► D^^D 11 channels 
shares the At background distributions with the com- 
binatorial background because it originates from similar 
sources. The B° — > D* + D*~ peaking background has 
only a lifetime component, since it originates from a spe- 
cific B + decay. Untagged events are also included in 
the fits to constrain the toes and cos#tr shapes but do 
not contribute to the determination of the CP parame- 
ters. We also allow the signal yield, the Toes background 
shape, and the cos^tr background shape to vary in the 
fits. Again we use high-statistics MC samples divided 
into data-sized subsets to validate the fitting procedure 
and find no significant bias. 

The statistical uncertainties of the CP measurements 
below are consistent with the expected uncertainties ob- 
tained from MC studies that include the signal and back- 
ground yields observed in data. The statistical uncer- 
tainty for the B° — > £)(*) ± £) T channels is essentially un- 
changed or even slightly worse than our previous mea- 
surement [11]. We interpret this as a downward fluctu- 
ation in the statistical uncertainty of the previous mea- 
surement. Using MC data, we estimate the probability 
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Af(ps) Af(ps) 

(c) B° -> D*+D~ (d) B° -* D+D*- 

FIG. 5: Projections onto At of the fit result and the data in the region rriES > 5.27GeV/c 2 for the three highest purity tagging 
categories. The triangular points and the dashed lines are for B° tagged events, and the circular points and solid lines are for 
B° tagged events. 



of observing such a fluctuation at about 20%. For each 
measurement that follows, the first uncertainty is statis- 
tical and the second is systematic. 

From the fit to the B° -> D* + D*~ data, we find 

S+ = -0.76 ±0.16 ±0.04 
C+ = +0.00 ±0.12 ±0.02 
S± = -1.80 ±0.70 ±0.16 
C_l = +0.41 ± 0.49 ± 0.08 , (11) 

with an effective R± — 0.155±0.030. If we perform the fit 
with the additional constraints that 5+ = S± = Sz>*+ rj»- 
and C+ = C± = C£>*+£>»-, we obtain 

S D ,+ D ,- = -0.70 ± 0.16 ± 0.03 

C D ,+ D - = +0.05 ± 0.09 ± 0.02 , (12) 



having an effective R± = 0.171 + 0.028. Fitting the B° -> 
D + D~ data yields 

S D+D - = -0.63 ± 0.36 ± 0.05 

C D+D - = -0.07 ± 0.23 ± 0.03 , (13) 

and fitting the B° —> D* ± Z? =F data yields 

S D , +D - = -0.62 ± 0.21 ± 0.03 
S D+D ,- = -0.73 ± 0.23 ± 0.05 
C D , +D - = +0.08 ± 0.17 ± 0.04 
C D+D *- = +0-00 ± 0.17 ± 0.03 . (14) 

Projections of the fit results onto At for events in the 
region tues > 5.27GeV/c 2 , and their flavor asymmetry, 
can be seen in Fig. 5. To enhance the visibility of the 
signal in these projections, we use three of the six tagging 
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categories with the highest purity, which account for 80% 
of the total effective tagging power Q. The correlations 
among the CP parameters are given in the appendix. 

We evaluate systematic uncertainties in the CP asym- 
metries for each mode by varying the fixed parame- 
ters for the mistag quantities and At resolution model 
within their uncertainties while accounting for correla- 
tions among the parameters. For the B° — > D+D~ and 
B° — > D*+D+ modes, we change the fixed ttt-es sig- 
nal width by ±0.2 MeV/c 2 , an amount determined from 
a comparison of data and MC event samples in modes 
with high purity, and take the difference in fitted results 
as a systematic uncertainty. Additionally, we vary the 
fraction and shape of the peaking background compo- 
nent. We also include systematics for possible detector 
misalignment and the presence of doubly-Cabibbo sup- 
pressed decays of the £?tag meson [32] . We assign a sys- 
tematic uncertainty equal to the statistical uncertainty of 
the MC sample used to validate the fit. Other sources of 
systematic uncertainty include: the B° meson properties 
(Am<j and t b o), which we vary to ±1(7 of their world av- 
erages, and uncertainty in the boost; the corresponding 
changes in the CP asymmetries are taken as the estimate 
of the systematic uncertainties. For the B° — > D*+D*~ 
mode, we vary the cos 8 tr resolution parameters and back- 
ground shape in the manner described for the evaluation 
of systematic uncertainties on R± and take the effects on 
the CP parameters as the associated systematic uncer- 
tainty. A summary of the systematic uncertainties for the 
CP parameters is given in Tables II and III. As before, 
the total systematic uncertainty is the sum in quadrature 
of the individual contributions. 

Because B° — ► D*+D+ decays are not CP eigenstates, 
it is illustrative to express the CP asymmetry parameters 
S and C in a slightly different paramctrization [33] 

Sd'D = ^ (Sd'+d- + S D + D ,-) 



ASd*d — g (Sd'+d- 



$d+D'- 



C 



D*D 



AC 



D*D 



2 {Cd'+d- 
2 {Cd'+d- 



C 



D+D* 



Cd+d*-) ■ 



(15) 



The Sd*d and Cd*d parameters characterize mixing- 
induced CP violation related to the angle (3 and flavor- 
dependent direct CP violation, respectively. ASd*d is 
insensitive to CP violation but is related to the strong 
phase difference 5. ACd»d describes the asymmetry be- 
tween the rates T(B° -> D*+D~) + T(B° -» D+D*-) 
and T(B° -> D+D*-) + T(B° -» D*+D~). Using the 
results from Eq. 14 and taking into account correlations 
among the variables, we find 

S D * D = -0.68 ±0.15 ±0.04 
ASd-d = +0.05 ± 0.15 ± 0.02 

C d *d = +0.04 ±0.12 ±0.03 
ACd*d = +0.04 ±0.12 ±0.03. (16) 



From the signal yields N D *+ D - and N D + D ,- deter- 
mined in the time-dependent fit described above, we also 
measure the time-integrated CP asymmetry in B° — > 
D* ± D+ decays, defined as 



A = 



N. 



D- + D- 



N D+D , 



N D , +D - + N D+D »- 



We find 



A = +0.008 ± 0.048(stat) ± 0.013(syst) . 



(17) 



(18) 



where the systematic uncertainty is dominated by track 
reconstruction efficiency differences for positive and neg- 
ative tracks (0.013). There is also a small contribution 
from the rriES signal width, peaking background, and MC 
statistics (0.002). 



V. CONCLUSION 

We have measured the CP asymmetry parameters for 
B° — > £)(*)+£)(*+ decays, including the CP-odd fraction 
in the B° -> D*+D*- channel, using the final BABAR 
data sample. All of the S parameters are consistent with 
the value of sin2/? measured in b ~ > (cc)s transitions [34] 
and with the expectation from the Standard Model for 
small penguin contributions. The C parameters are con- 
sistent with zero in all modes. In particular, we see no 
evidence of the large direct CP violation reported by the 
Belle Collaboration in the B° — ► D+D- channel [15]. 
This measurement supersedes the previous BaBAR mea- 
surements [11, 12] of CP asymmetries in these decays. 
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TABLE II: Systematic uncertainties on the B° — > D* + D* CP paramefc 


ers. 






S+ 


S± 


c+ 


Cx 


V + D«- 


C D *+ D *- 


Tagging and At resolution 


0.022 


0.031 


0.010 


0.017 


0.021 


0.009 


Peaking background 


0.012 


0.079 


0.002 


0.019 


0.012 


0.003 


Detector Alignment 


0.006 


0.029 


0.001 


0.019 


0.005 


0.002 


Doubly-Cabibbo suppressed decays 


0.002 


0.002 


0.014 


0.014 


0.002 


0.014 


Potential Fit Bias 


0.011 


0.098 


0.008 


0.065 


0.011 


0.007 


Angular PDF variations 


0.025 


0.091 


0.004 


0.015 


0.011 


0.001 


Other 


0.013 


0.025 


0.005 


0.029 


0.013 


0.002 


Total 


0.040 


0.163 


0.020 


0.080 


0.032 


0.018 



TABLE III: Systematic uncertainties on the B° -> D M± D* CP parameters. 





$d+d- 


C D + D - 


Sd*+d- 


C>D* + D- 


Sd+d*- 


Cd+d*- 


Tagging and At resolution 


0.031 


0.011 


0.027 


0.012 


0.029 


0.011 


toes signal width 


0.034 


0.020 


0.013 


0.018 


0.028 


0.012 


Peaking background 


0.018 


0.007 


0.014 


0.023 


0.030 


0.013 


Detector Alignment 


0.002 


0.001 


0.004 


0.002 


0.002 


0.001 


Doubly-Cabibbo suppressed decays 


0.002 


0.014 


0.002 


0.014 


0.002 


0.014 


Potential Fit Bias 


0.007 


0.005 


0.008 


0.006 


0.008 


0.006 


Other 


0.006 


0.002 


0.002 


0.006 


0.004 


0.003 


Total 


0.051 


0.028 


0.034 


0.036 


0.051 


0.026 



TABLE IV: Correlations among the CP parameters of the 
B° -> D* + D*~ mode split by CP-even and CP-odd. 





s+ 


c+ 


s ± 


Cx 


Px 


s+ 


1 


0.008 


0.376 


-0.036 


-0.083 


c+ 




1 


0.045 


-0.465 


0.003 


s± 






1 


-0.224 


0.471 


Cx 








1 


-0.151 


R± 










1 



APPENDIX: CORRELATIONS AMONG THE CP 
PARAMETERS 

To allow detailed use of these results, we include the 
correlation matrices for the CP parameters. Table IV 



contains correlations among the fit parameters in the 
B° — > D* + D*~ channel with separate CP-even and CP- 
odd asymmetries, and in the combined case, the corre- 
lation between S D *+ D *- and C D »+ D »- is 0.8% with cor- 
relations to the effective Px the same as the CP-even 
parameters. Table V contains the correlations among 
the B° — > D*^D T asymmetries. The correlation of the 
time-integrated CP asymmetry A with any of the CP 
parameters is less than 0.1%. The correlation between 
Sd+d- and C D + D - is —1.2%. 
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